The internal states of active inclusions and the dynamics of an active membrane 
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A theoretical model of a two-component fluid membrane containing lipids and two-state active 
inclusions is presented. Under strong inclusion activities the membrane can be unstable due to pump- 
driven undulation or aggregation of high curvature excited inclusions. Depending on the structural 
details of the inclusions, the active conformation change of the inclusions can also induce finite 
size domains on the membrane. At long wavelengths a stable membrane has a height fluctuation 
spectrum similar to previous studies which neglected the inclusion internal states. For a tensionless 
membrane, (h(q ± ,t)h(—q±,t)) ~ TefjqJ 4 , where T e // is an effective temperature which depends 
on the details of inclusion activities. 

PACS numbers: 87.16.-b, 05.40.-a, 05.70.Np 



A biological membranes is a multi-component mixture 
of lipids, proteins, carbohydrates, and other materials. Q 
Since the active inclusions (e.g., proteins) in the mem- 
branes participate in physiological processes such as ion 
transport, signal transduction, and cell locomotion, re- 
cent theoretical and experimental studies have focused 
on the collective effect of these active inclusions on the 
membranes. 0, 0,13 An important conclusion of these 
studies is that, unlike equilibrium membranes, [f| the 
height fluctuation of a tensionless active membrane in the 
long wavelength regime depends on the inclusion activ- 
ities. Activity-induced instabilities are also predicted in 
previous theoretical models. 0,0 On the other hand, the 
fact that these active inclusions have more than one in- 
ternal conformation state has been discussed for a single 
transmembrane protein and for a single ion pump || . 
Although the elastic interaction between lipids and an 
active inclusion have significant effect on the functions of 
a single inclusion, the equally important question of how 
this interaction affects the collective behavior of the lipid- 
inclusion system is seldom addressed in previous studies 
of active membranes. 

In this Letter, I discuss the dynamics of a fluid mem- 
brane containing lipids and identical two-state active in- 
clusions. As shown in Figure 1, inclusions in different 
conformation states have different couplings to the local 
mean curvature of the membrane. Recent experiments on 
the structures of potassium channels 01 an d K + chan- 
nels have provided structural evidences of such inclu- 
sion conformation changes. A conformation change of an 
inclusion can be induced by external stimuli (active tran- 
sition) or by the lateral pressure exerted by the surround- 
ing medium (passive transition). During a conformation 
change, an inclusion also exerts forces to the lipids and 
the solvent. By taking these effects into account, the 
main predictions of the model are: (i) The height vari- 
ance of a stable tensionless membrane in the long wave- 
length limit is (\h(q±, t)\ 2 ) ~ T e f / / 'q\, where T e ff is an 



effective temperature which depends on the inclusion ac- 
tivities. For a stable tense membrane, surface tension 
dominates the long wavelength membrane height vari- 
ance, therefore (\h(q±,t)\ 2 ) ~ T jq\_, the same as that of 
an equilibrium tense membrane, (ii) The system has two 
types of long wavelength instabilities. They are pump- 
driven undulation instability due to the pumping-induced 
attraction between the inclusions, and the excitation- 
driven undulation instability due to excitation-induced 
aggregation of the excited inclusions in the case when 
the excited inclusions have sufficiently higher curvature 
than the ground state inclusions, (iii) Depending on the 
structural details of the inclusions, the system can have a 
finite wavelength instability which leads to domains with 
a characteristic length depending on the active transition 
rates. At high active transition rates, the activities can 
even drive the system to a homogeneous state. There- 
fore it is possible to control the membrane morphology 
in experiments by tuning the activities of the inclusions. 

In the Monge representation the shape of the mem- 
brane is described by h(x,y). The number density of 
inclusions in state a (a = 1 or 2) is <fr a (r_L,i) where 
rj_ = (x,y). It is useful to introduce $± = (<f>i ± <& 2 )/2 
and = J d 2 r±$ + / J d 2 r . The coarse grained Hamil- 
tonian of the system H = H m + Hi + H c includes the elas- 
tic energy of the membrane, the direct interaction energy 
of the inclusions, and the membrane-inclusion coupling. 
To lowest order H m has the form 

H m = i J d 2 r ± (K(Wlh) 2 + j(W ± h) 2 ) , (1) 

where n is the bending rigidity, 7 is the surface tension 
of the membrane. In this Letter I consider the regime 
where the system is close to a phase separation in 
and <&+ is non-critical. Therefore the interaction energy 
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of the inclusions in the harmonic theory is 0, 0] 
Hi = J d 2 r ± { I(£ 2 (V ± $_) 2 + r 2 $ 2 _) 



in 



2 (* + _$fy»_A 6 *_}, (2) 

where r 2 > is small and to ;g> r 2 . A e > is the exci- 
tation energy for an inclusion to change its conformation 
from state 1 to state 2. Thus "state 1 (2)" is the "ground 
(excited) state" of an inclusion. The couplings between 
the inclusion densities and the membrane curvature de- 
pend on the shapes of the inclusions, it can be expressed 
by 

H c = J d 2 r± _ k(ci$i + c 2 $ 2 )V 2 L h 

= J d2r ±- «(c+$+ + c_$_)Vift, (3) 

where c± = c\ ± c 2 . The equilibrium properties of 
this model Hamiltonian is briefly summarized in the 
following. The stability criterion of the membrane is 

2 > 0. For a ten- 
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ca(g±) = («- - j^+ r2i 

sionless membrane there is a long wavelength instability 
when 



K eff = « - 
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<0, 
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where a phase separation in $_ occurs. For a membrane 
under tension the system has a finite wavelength insta- 
bility when 



< 



7— + 2a K+ (j— 

, r 2 V V f2; 
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where k + 



cr\/m. In this case, the long wave- 



length instability is cut off by the surface tension at 
large length scales. When the system is homogeneous, 
the equilibrium height fluctuation of the membrane in 

for a 
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for 



the long wave limit is < h(qj_)h(—qj_) >= ; ' 

tensionless membrane and < h(q±)h(—q±) >= 
a membrane under tension. 

To discuss the dynamics, I introduce suitable variables 
4>±. When c 2 {q±) > 0, </>± are chosen to be the devia- 



tion from the equilibrium solution, i.e., 



(0) 



and 



(0) 



value of $± 
0+ = $+ 



where $^ 



is the equilibrium 



$^ 0) and , 



When c 2 {q±) < 0, 0± are chosen to be 
= such that the effect of non- 
equilibrium activities on the stability of the membrane 
can be discussed conveniently. On most experimental 
time scales the membrane does not exchange inclusions 
with the solvent, therefore <fi + obeys the conserved dy- 
namics 



at 



A. 



" ± d<t> + 
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where the vector thermal noise C+ has zero mean 
and variance (Vj_ ■ C + (r,t)V ± ' C+( r \t')) = 
-2fc B TA + V 2 <5 3 (r - r')S(t - t') . In general 0_ 
satisfies an equation of the following form, 

r)A> „ SH 

= (-A + A 2 Vi)— -fe + 0+-fe_^_+C- (7) 



at 



The passive conformation change of the inclusions is 
taken into account by Ao term, 14] A 2 term represents 
the mutual diffusion of the inclusions, the conformation 
changes due to external stimulations (active transitions) 
are taken into account by — k + 4> + — k—4>- the ther- 
mal noise has zero mean and variance 2ksT(\Q + 
^2lj_)- The hydrodynamic flow, the permeation of sol- 
vent through the membrane, and the force exerted on the 
membrane by the inclusions during active transitions all 
contribute to the dynamics of the membrane, thus 



dh , (SH 



+ p;m>2 )+Ch, (8) 



here P£ (P^) is the momentum transferred from an inclu- 
sion during an active excitation (relaxation) process. The 
thermal noise ^ has zero mean and variance 2ksTX p . 
The solvent which embeds the membrane satisfies the 
modified Stokes equation 

SH 

= -Vp(r,t)-— (v ± ,t)S(z)z 

+P°[S(z - <) - 5{z + wf))k e (f)iz 

+ Pa[S(z - Wl) - S(Z + wi)]k r fci 
+7 ? V 2 V + f 1 ,. 



(9) 



Here p is the pressure, 77 is the solvent viscosity, the third 
and forth terms on the right hand side come from the 
active transitions of the inclusions. Although the total 
force acting on the solvent due to an inclusion confor- 
mation change is zero, the force distribution is not van- 
ishing. To lowest order we approximate the force dis- 
tribution of an active inclusion excitation (relaxation) 
event by a force dipole. |3| w", wf, w 2l and w 2 are 
characteristic lengths associated with inclusion confor- 
mation transition processes. The vector thermal noise 
ft, has zero mean and variance (f V i(r,t)f V j(T',t')) = 
2k B TT 1 {-S vj V 2 + didjJS 3 ^ - r')S(t ~ t'). 

Solving the Stokes equation @ and consider the 
regime where the dynamics of 0+ is fast compared to 
that of 4>-, 16], the linearized equations of motion in the 
Fourier space for h and </>_ are obtained. Since the contri- 
bution from permeation in a typical system is negligible 
for lengths small compare to O(lcm) Q, the equations 
for h and <f>_ becomes 



dh(g ± ,t) 

04>-(q±,t) 
at 



D h (qx) Dh<t>{qx) 

D$h(qx) P>$(qx) 

+ ( fh(q±,t) 

V U (?-!>*) 



h(qx,t) 

<f>-{qx,t) 



(10) 
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where 



1 



4f)q± \\ m 
D4q±) = -(A + A 2 gi)(^ 2 gl + r 2 ) - fc. 

Dh<p(q±) 
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and 



D<f,h(q±) = (Ao + A 2 gi)KC_(7j 



fc + ^±,i, (11) 
m 

kcI ee kc± + u±, and u± = P Q e fc e (-(u^) 2 + «) 2 )/2 ± 
P£fc r (-(w 2 ) 2 + (wi) 2 )/2. It is convenient to think c% as 
the "renormalized" inclusion-membrane coupling in the 
presence of non-equilibrium activities. Magnitudes of v± 
depend on structural details and active transition rates 
of the inclusions. Dimensional analysis in previous stud- 
ies 0,0 show that v± ~ fl 2 , where I ~ 5nm is a typical 
length associated with an inclusion, and / has the di- 
mension of a force. At high activities / ~ 10~ 12 N 0,0, 
together with nc± ~ IksT gives v±/kc± < 0(1). There- 
fore kcI and kc+ (kc" and kc_) are of the same order of 
magnitude, and usually have the same sign. The noises 
fh(q±,t) and f$(q_\_,t) have zero mean and variances 
(A(<U,*)/fc(q'x,f)} = ^(l + v 2 + q ± /ir]X + m 2 )6(t- 
t')5 2 (q±_ + q'J, and (^(q x , t)/*(q' ± , i')) = 
2fc B T(A + X 2 ql + k 2 + /\+m 2 ql) S(t - t')S 2 (qx + q'j.)- 
Notice that there are contributions from the activities of 
the inclusions. They are no longer equilibrium thermal 
noises. 

The stability of the system is determined by the eigen- 
values of the 2x2 matrix in Eq. lfTT))l . When both Dh and 
Dff, are positive, the stability condition for a tensionless 
membrane is 



+ [Xo< ff + 



k + k 2 c + c°l 
r 2 m 



—k1 > 0. 

r 2 J 



(12) 



where k° = k — K 2 c + c a + jm, K®ff = 
It immediately follows that the system has a long wave- 



length instability when v+- 



is positive and 



sufficiently large such that the re"^ is large and neg- 
ative. This instability is due to the pumping of the 
active inclusions. As shown in Figure 2, when active 
pumping attracts more inclusions, becomes nega- 

tive. A similar instability is discussed in a previous study 
which neglected the internal states of the inclusions [2j. 
Another type of long wavelength instability occurs at 
n 2 c + ct < and sufficiently large k+ = k e — k r . Since ct 
is the renormalized c_ in the non-equilibrium situation, 
k 2 c+c1 ~ k 2 c+c_ < can be interpreted intuitively as 
the case when c 2 2 > ci 2 , i.e., the excited state prefers 
higher membrane curvature. Thus, as Figure 3 shows, 
this instability is due to the strong active excitations of 
the inclusions to their higher curvature state. 



The system has a finite wavelength instability when 
the q\ and q 1 ^ terms in Eq. I|12|) are positive and 
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/, k+ K 2 c+c a fc_ 

I A 2 K" X K a eff + — — + — K% 

/r 2 U + V r 2 m r 2 + 



(13) 



The resulting steady state has an average domain size 
R which can be estimated by the following scaling anal- 
ysis. 0, When there is no budding, the growth of 
intramembrane domains in the absence of inclusion ac- 
tivities corresponds to a two-dimensional phase separa- 
tion dynamics, i.e., R(t) ~ t a , a sa 1/3 [H- This 
growth eventually saturates due to the active transi- 
tions. Since the time scales associated with active tran- 
sitions are k±, the steady state domain size R should 
obey R ~ k^f, where fcju = when k+K 2 c + ct/m is 

greater (smaller) than k-K°L. Notice that this instabil- 
ity is suppressed at large fc_ = k e + k r . This is because 
when the inclusions change their conformations at ex- 
tremely high rates, the membrane feels the time-averaged 
inclusion conformation, therefore there is no finite size 
domains in the steady state. When k 2 c + c°l > 0, an ex- 
cited state inclusion has smaller curvature than a ground 
state inclusion, sufficiently large fc + stabilizes the mem- 
brane because strong active transitions excite the inclu- 
sions to the smaller curvature excited state. 

When the membrane is stable, its height fluctuations in 
the long wavelength limit is expressed by ^ *ff , where 



T, 



Tk 



eff 



eff 



(14) 



r 2 -fc-/A 



If the passive conformation changes of the inclusions are 
extremely rare such that Ao <C X 2 q\ in all experimental 
length scales, T e // has the following form, 



T eff ~ 



Tk 



eff 



(15) 



Eqs. I|14 |l l|15 |l are similar to previous studies on the active 
membranes which neglected the internal states of the in- 
clusions. Ijl For a membrane under tension the kBT \ ss 

behavior is cut off at long wavelengths. The membrane 
height fluctuation at large lengths in this case is simply 
Y7<Zj_ where T is the temperature of the solvent. 
In summary, a theory of active membranes with two- 
state active inclusions is discussed. Although current 
analysis is restricted to a relatively small region of the 
parameter space, the rich dynamical behaviors shown in 
the analysis are believed to exist in more general situa- 
tions. These include instability induced by active pump- 
ing, instability induced by inclusion-excitation, finite size 
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domains controlled by activities, and activity-restored 
homogeneous morphology. In the stable case, the mem- 
brane height fluctuation is similar to several previous the- 
oretical and experimental studies. The analysis of this 
model in a broader parameter range and the effects of 
nonlinearity will be included in a future work. |2l| It is 
my hope that this model might be of relevance to ac- 
tive conformational transitions in biological membranes, 
and experimental work on artificial membranes might be 
designed to confirm the main ideas of this model. 
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FIG. 1: Schematics of an inclusion in different conformation 
states 
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FIG. 2: Large negative n^ff instability: pumping-induced 
attraction between inclusions induces inclusion aggregation. 
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FIG. 3: The presence of many high-curvature excited state 
inclusions induces a long wavelength instability. 



